for i in [2, 2, 4, 8, 16, 32, 64]
    T1 = time()
    global nels = i
    include("main1.jl")
    println(neq, "   cpu time = ", time() - T1)
end
